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ABSTRACT 


The problem of trajectory optimization is considered from the stand- 
point of numerical analysis. A numerical solution is obtained for an 
assumed thrust angle history, and an explicit numerical solution is 
obtained for the linearized equations of motion for neighboring solutions* 
With the explicit solution available, it is not difficult to determine 
whether or not the assumed solution is optimum. The first and second 
variations together provide a straightforward iteration method of approach- 
ing the optimum solution. It is of particular interest that the procedure 
remains unaltered even by the introduction of discontinuities and inter- 
mediate constraints . 
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DEFINITION OF SYMBOLS 


Symbol 


Definition 


X = 



X 


fa 
gx 
g2 
hi 
h2 
m x 
m 2 J 


> 


L n 


t 


time derivative of X 


impulse response functions defined by equations (4) and (12) 


radial distance 

required radial distance at cutoff time 
actual radial distance at cutoff time 
time 

cutoff time 
velocity 

velocity at cutoff time 


Cartesian coordinates of position in a space fixed 
system 


iii 


Symbol 

• • 

x,y 


z 


Zi 


Ar 

At 

Ay 


AA^, AA^j AA3 
AX. 


\l» ^2> ^3 
\l> ^2j ^3 

X 

X 


DEFINITION OF SYMBOLS (Continued) 

Definition 

time derivatives of x and y, respectively 

an arbitrary function of X and t 

the value of z at time t x 

required value of z at t x 

deviation in cutoff radial distance 

deviation in cutoff time 

deviation in cutoff velocity 

corrections for A x , A 2 and A3, respectively 
variation in thrust angle 
correction to X 

coefficients determined such that h x + A x f x + A^gx + A^mx = 0 
first approximations to A x , A 2 and A3, respectively 
thrust angle measured from the y-axis 
first approximation for X 


iv 
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TRAJECTORY OPTIMIZATION BY EXPLICIT NUMERICAL METHODS 

SUMMARY 


The problem of trajectory optimization is considered from the stand- 
point of numerical analysis. A numerical solution is obtained for an 
assumed thrust angle history, and an explicit numerical solution is 
obtained for the linearized equations of motion for neighboring solutions. 
With the explicit solution available, it is not difficult to determine 
whether or not the assumed solution is optimum. The first and second 
variations together provide a straightforward iteration method of approach- 
ing the optimum solution. It is of particular interest that the procedure 
remains unaltered even by the introduction of discontinuities and inter- 
mediate cons traints . 


OPTIMIZATION WITH END POINT CONSTRAINTS 


The problem of trajectory optimization can be considerably simplified 
in concept if advantage is taken of available numerical methods from the 
beginning. To illustrate the point, the problem of minimizing the burning 
time required to meet desired end conditions will be considered. The dif- 
ferential equations describing the motion of the vehicle are defined in 
the following vector equation: 


X = F(X,X, t). 


( 1 ) 


The problem is to minimize the time required to fulfill the following 
conditions at cutoff: 


r (t c ) = r c 

e(t c ) = e c 

v(t c ) = v c 


( 2 ) 


For some function X(t) with a velocity cutoff constraint, the follow- 
ing numerical solution is obtained: 



(3) 


The linearized equations of motion can be solved explicitly for the 
same initial conditions and cutoff conditions to give a solution of the 
following form*: 




^ = r n “ r c + 


fj/CX + faAX 2 + .. 


dt 


^n, 


A9 = e n - e c + 


gl AK + g 2 A)C 2 + 


dt 


At 


L. 

- r 


n r- 


J 

t 


hjAX + h 2 AX 2 + ... 


dt 




(4) 


where 



At - t - t n 
AX - X(t) - X(t) 


J 


The derivation of these equations is described in detail 


in Reference 2 . 
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The above equations can be combined in the following expression: 


At + A x Ar + ApAB = 



(hx + Axf i + Apg i ) AC 


+ 


(h-2 + Axf 2 + Apgp) AX. 2 + 


dt + Ax(r n - r c ) 


+ A 2(e n - e c ). 


This expression is general at this point and imposes no constraint on 
r, 0 or AX.(t). Setting r = r c and 0 - 0 C gives the following expression 
for At: 


At 


L 

■/ 


n 


(h x + A x f x + Apg i ) AX + (h 2 + A x f 2 + Apgp) 


AX 2 + 


dt 


+ A x (r n - r c ) + A 2 (6 n - 0 C ). 


( 6 ) 


If h x + A x f x + Asg x = 0 and (h 2 + A x f 2 + Aag 2 ) > 0, any change in the 
function AX(t) over any interval will increase the value of At in a 
neighborhood. If this condition exists, the function X(t) yields a 
minimum value of t locally. If h x + A x f x + Apg n ^ 0, consider A x and A 2 
such that 


h x + A x f x + A 2 g x - e. 


The values of A x and A 2 can be obtained by the method of least squares or 
any other comparable method. Setting A x = A x + AA X , A 2 = A 2 + AA 2 and 
AX = A£ + SX gives the following expression for At. 
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At = 


e (t) + 


AA x fx + 


AAgg i + 2(1*2 "*■ ■*" Apgp) AX 


SX + 


-} 


dt 


+ ^i( r n ” r c ) + ^( e n " e c ^* 


,(7) 


Setting the weighting function for SX equal to zero gives the 
following equation: 


2 (h 2 + A^fg + Apg p) AX — -e(t) - fxAAx “ 


( 8 ) 


For simplification, the following definitions are employed: 


h 2 - (h 2 + A x f 2 + A2g 2 ) A 


e(t) = + 


ejtl 

2h| 


fi - 


2h* 


" 2h* 


(9) 


J 


The expression for AX is the following: 


AX — f xAAx i g xAA 2 — e . 


( 10 ) 
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The terms AA X and AA 2 can be determined by setting Ar = 0 and A0 = 0 in 
equations (4) and keeping only linear terms. This results in the follow- 
ing system of equations: 



( 11 ) 


Solving this system for AA X and AX 2 provides the necessary quantities 
to determine AX from equation (10) at each time point along the trajectory. 
With the new function X x = X + AX, another iteration can be performed to 
determine X 2 . This can be continued until r n - r c , 9 n - 9c and e(t) all 
meet the required tolerances. Convergence of this iteration process 
ensures that the first order weighting function in equation (6) is zero. 
Similar operations can be applied to the second order term in the event 
that it is not everywhere positive. 


INTERMEDIATE CONSTRAINTS 


Assume that at some other time, t x , which may or may not be a func- 
tion of X, an additional constraint is imposed requiring that some other 
function z(t x ) = z c and that for X(t), z(t x ) = z x . The expression for Az 
can be written as 


Az = z x - 


m x AX + H 12 AX 2 + 


dt. 


( 12 ) 


If t x > t n , m x = 0 and m 2 = 0 for tn < t § ti. Equation (12) can 
then be written as 


Az = z x - z c + 



m x AX + m^AX 2 + 


dt. 


(13) 
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If to < t x S t n , define mi = m 2 = 0 for t x < t § t n . Then Az can 
also be represented by equation (13). This equation is identical in 
form to expressions for Ar and A 0 appearing in equations (4). The same 
analysis used before can be applied to the constraints that 


r(t c ) = r c 
9(t c ) = 9c 
v(t c ) = v c 
z(ti) = z c . 


The resulting expression is obtained for Z^X: 


f § 1^2 TO2AA3 4" e 


(14) 


where 


fi - - — ~ , h’^ - h 2 + A x f 2 + Xgg 2 + 


gx - 


®i = 


2h. 


2h* 


m-, 


2h 


Vc 


€ = 


2h: 


(15) 


and 7\j_, A 2 and A3 have been determined as before. The values AA x , AA2 
and AA3 are solutions to the following system of equations: 


>n 


rn 


AA 


1 / f i f i dt + AA 2 / fig! dt + AA 3 / fifii dt = / f x e dt + r c - r n 


to 
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AAn 


tn t 

/ gJidt + AA2 I Sail dt + M 3 J 


n 


gi®]. dt = 


gl e dt + e c - e n 


AAi 


t n t n t n t n 

J dt + AA 2 J m iSi dt + AA 3 ^ m-jfi-L dt = J' m x e dt + z Q - z x 


The fact that the form of this system remains the same whether the con- 
straints are applied to the end conditions or at some other point along the 
trajectory makes this method readily adaptable to computer operations and 
an extremely versatile multistage optimization program could easily be 
developed with a wide choice of constraints imposed at arbitrary points 
along the trajectory. 



APPENDIX 


An example to illustrate the method has been completed* from the 
following initial conditions: 


: = 236.25546 km 
o 

x Q = 2800.5746 m/sec 

■ 0 = 6456.66782 km 

y Q = 635.94317 m/sec 

: Q = 184.99481 sec. 



The motion is assumed to be described by the following two- 
dimensional system of differential equations: 


x 


F . v , 

— sm X 4- x , 
m g 


y 


F 

m 


cos X + y , 



y 



g> 


g 


_ 2 

g r 
o o 


g Q = 9.81 m/sec 2 , 


r = 6370 km 
o 


F _ 10,787315 = t - 184.99481 

m 1.2867664- . 261282644 t x ’ Tl 100 


184.99481 S t < 412.85503 


F = 8.92405155 = t - 412.85503 

m .69140715 - . 212616813 t 2 ’ 2 100 

412.85503 g t < 523.87973 


- = 0, 523.87973 g t < 528.87973. 

m 


The computer program used to solve this example was developed in coopera- 
tion with James Hilliard, R-C0MP, MSFC, and is described in Reference 4. 
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preceding page 


blank, not filmed. 


F _ 1.96132994 t - 428.87973 

m .36021199 - .046948359 t 3 » T3 " 100 

t £ 428.87973. 


F/m is in m/sec 2 and t is in seconds. 

The problem is to determine the function X which will minimize 
the burning time required to reach the following conditions: 

r = r = 6555.2 km 
c 

v = v c = 7792.5746 m/sec 

0=0 = 90°. 

c 

In addition to the above constraints, an intermediate constraint 
has been included which requires that the spent stage which separates 
at t = 412.85503 seconds must, moving only under the influence of 
gravity, pass through the prespecified point (x x , y x ) where 


Xi = 2588.9871 km 
y z = 5971.13760 km. 


The following function was initially assumed for X: 


X = 75.111 + 1.1982 


where X is in degrees and t in seconds. 

The equations of motion were solved numerically by fourth order 
Runge Kutta at five-second time steps until the stage dropping off at 
t = 412.85503 seconds fulfilled the condition that x = x z and the 
upper stage fulfilled the condition that v = v c . The weighting 
functions, f z , g x , m x and k x , were computed numerically by the method 
outlined in Chapter I of Reference 2. Corrections to the function X 
were determined from equations (14) and (15), and another iteration 
was begun. After the seventh iteration, the following conditions were 
fulfilled: 
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r(t c ) 

v(t c ) 

e(t c ) 


r c ± 0.1m 
v c ± 0.01 m/sec 
0 C ± 10“ 5 radians 


xCtx) = x x ± 0.1m 

y(t x ) = Yi ± o.lm 


e(t^) = 0 ± 10“ 6 / radian. 


Figure 1 shows the resulting trajectory compared with the solution 
that resulted when the intermediate constraint was removed. The total 
burning time with the additional constraint required 768.60 seconds com- 
pared with 690.24 seconds without this constraint. The additional 78.36 
seconds indicates that this is a rather expensive maneuver but illustrates 
the versatility of the method. 

The guidance angle, X, is shown in Figure 2, where it is compared with 
the function used for the first iteration. A complete tabulation of X at 
every point for which it was computed is included in Table 1 at the end 
of the appendix. In performing the numerical integration, this table 
was used with second order interpolation employed to determine the mid- 
point value required by the Runge-Kutta method. 

Figure 3 illustrates the first order weighting functions f x , g-L, m ls 
and f 1# These functions are tabulated for several time points in Table 1. 
For simplicity they have been rounded off to the fourth decimal place. 

It can be readily verified that they satisfy the equation h x + 'hjfx + 

Aag x + A^mj = 0 within a tolerance of less than one unit in the last 
significant figure for the following values of A x , A 2 and A 3 . 


7\ x = .299026 sec/km 

A 2 = .702527 10 3 sec/radian 

A 3 = -.633224 sec/km. 


Figure 4 shows the functions f l5 g x and h-L which were used to deter- 
mine the corrections to X. The second order weighting function, h 2 , is 
also shown in this illustration, and it can be seen to be positive. 
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FIG. 2. v FOR THE INTERMEDIATE CONSTRAINT TRAJECTORY 
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TABLE 1 


t(sec) 

f ifkm/ rad 

sec] g 1 [10“ 3 /sec] 

Hi! [km/ rad sec] 

h^l/rad] 

185 

-2.2642 

.9238 

-1.9664 

-1.2171 

200 

-2.2895 

.9679 

-1.9963 

-1.2594 

250 

-2.3827 

1.1380 

-2.1054 

-1.4201 

300 

- 2.4932 

1.3542 

-2.2331 

-1.6199 

350 

-2.6290 

1.6382 

-2.3880 

-1.8769 

400 

-2.8032 

2.0267 

-2.5837 

-2.2217 

412.9 

-2.8568 

2.1507 

-2.6433 

-2.3305 

412.9 

-4.8807 

1.8953 

0 

.1280 

450 

-4.9237 

2.1708 

0 

- .0527 

500 

-4.9454 

2.6950 

0 

- .4146 

523.9 

-4.9455 

3.0395 

0 

- .6565 

523.9 

0 

0 

0 

0 

528.9 

0 

0 

0 

0 

528.9 

-1.3418 

.8472 

0 

- .1939 

550 

-1.2401 

.8836 

0 

- .2499 

600 

- .9809 

.9823 

0 

- .3968 

650 

- .7032 

1.1000 

0 

- .5625 

700 

- .4169 

1.2384 

0 

- .7454 

750 

- .1317 

1.3997 

0 

- .9440 

773.6 

.0000 

1.4886 

0 

-1.0458 


> 

H 

II 

.299026 sec/km 




a 2 - 

.702527 10 3 sec/rad 




A 3 = -.633224 sec/km 
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TABLE 2 


t(sec) 

X(degrees) 

t (sec) 

X(degrees) 

184.99481 

50.0558152 

325 

54.7621191 

185 

50.0559990 

330 

54.9239170 

190 

50.2312133 

335 

55.0854054 

195 

50.4058228 

340 

55.2465994 

200 

50.5798353 

345 

55.4075148 

205 

50.7532584 

350 

55.5681676 

210 

50.9261002 

355 

55.7285747 

215 

51.0983688 

360 

55.8887534 

220 

51.2700725 

365 

56.0487214 

225 

51.4412196 

370 

56.2084973 

230 

51.6118189 

375 

56.3681001 

235 

51.7818790 

380 

56.5275498 

240 

51.9514090 

385 

56.6868666 

245 

52.1204181 

390 

56.8460719 

250 

52.2889156 

395 

57.0051875 

255 

52.4569111 

400 

57.1642363 

260 

52.6244145 

405 

57.3232420 

265 

52.7914357 

410 

57.4822289 

270 

52.9579851 

412.855030 

57.5733702 

275 

53.1240733 

412.855030 

125.234677 

280 

53.2897111 

415 

124.935432 

285 

53.4549096 

420 

124.236574 

290 

53.6196802 

425 

123.535040 

295 

53.7840346 

430 

122.830917 

300 

53.9479848 

435 

122.124289 

305 

54.1115432 

440 

121.415244 

310 

54.2747226 

445 

120.703869 

315 

54.4375359 

450 

119.990252 

320 

54.5999968 

455 

119.274482 
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TABLE 2 (Continued) 


t(sec) 

X(degrees) 

t(sec) 

X(degrees) 

460 

118.556646 

585 

100.202874 

465 

117.836835 

590 

99.4630176 

470 

117.115138 

595 

98.7235515 

475 

116.391644 

600 

97.9845716 

480 

115.666444 

605 

97.2461745 

485 

114.939627 

610 

96.5084576 

490 

114.211284 

615 

95.7715192 

495 

113.481505 

620 

95.0354585 

500 

112.750381 

625 

94.3003754 

505 

112.018000 

630 

93.5663709 

510 

111.284453 

635 

92.8335468 

515 

110.549831 

640 

92.1020062 

520 

109.814223 

645 

91.3718527 

523.879730 

109.242771 

650 

90.6431911 

523.879730 

109.242771 

655 

89.9161274 

526.379730 

108.874254 

660 

89.1907682 

528.879730 

108.505553 

665 

88.4672214 

528.879730 

108.505553 

670 

87.7455956 

530 

108,340267 

675 

87.0260004 

535 

107.602271 

680 

86.3085466 

540 

106.863648 

685 

85.5933453 

545 

106.124489 

690 

84.8805090 

550 

105.384885 

695 

84.1701506 

555 

104.644926 

700 

83.4623837 

560 

103.904703 

705 

82.7573229 

565 

103.164310 

710 

82.0550831 

570 

102.423837 

715 

81.3557796 

575 

101.683378 

720 

80.6595283 

580 

100.943026 
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TABLE 2 (concluded) 


t(sec) 

725 

730 

735 

740 

745 

750 

755 

760 

765 

770 

773.599673 


X(degrees) 

79.9664454 

79.2766473 

78.5902504 

77.9073712 

77.2281261 

76.5526312 

75.8810022 

75.2133544 

74.5498023 

73.8904599 

73.4188071 



REFERENCES 


1. Bellman, Richard, "Introduction to Matrix Analysis," Chapter 10, 
McGraw-Hill, New York, 1960. 

2. Dickey, Lyle R. , "Guidance Applications of Linear Analysis," 

NASA TM X-53166, MSFC, November 27, 1964. 

3. Greenspan, Donald, "On Approximating Extremals for Boundary Value 
Problems," MRC Technical Summary Report #466, The University of 
Wisconsin, Madison Wisconsin, March 1964. 

4. Hilliard, James, "Computer Program: Trajectory Optimization by 

Explicit Numerical Methods," to be published as a NASA TM X. 


20 


APPROVAL 


NASA TM X- 53464 


TRAJECTORY OPTIMIZATION BY EXPLICIT NUMERICAL METHODS 
By Lyle R. Dickey 

The information in this report has been reviewed for security 
classification. Review of any information concerning Department of 
Defense or Atomic Energy Commission programs has been made by the MSFC 
Security Classification Officer. This report, in its entirety, has 
been determined to be unclassified. 

This document has also been reviewed and approved for technical 
accuracy. 


/S/ -* . ^ 

E. D. Geissler 

Director, Aero-Astrodynamics Laboratory 


21 


DISTRIBUTION 


R-DIR 

DEP-T 

I-I-IB-DIR 
Col. James 

I-MO-MGR 
Dr. Speer 

I-MO-O 
Mr. Kurtz 

R-RP-U 

Roy Bland 

MS- IP 

MS-IPL (8) 

MS-H 

CC-P 

MS-T (6) 

R-P&VE 

Mr. Swanson, -J 
Mr. Goerner, -A 

R-ASTR 

Mr. Moore 
Mr. Richard 
Mr. Gassaway 
Mr. Taylor 
Mr . Brooks 
Mr. Hosenthien 
Mr. Scofield 
Mr. Woods 
Mr. Digesu 
Mr. R. Hill 
Mrs. Neighbors 
Dr. R. Decher 


R-AERO 

Dr. Geissler 
Mr. Jean 
Mr . Dear-man 
Mr. Schwaniger 
Mr . Ingram 
Dr. H. Krause 
Mr. R’ne infur th 
Mr . Thomae 
Mr . Baker 
Mr . Burns 
Mr. W. Vaughan 
Mr. Hart 
Mr. Lovingood 
Mrs . Chandler 
Mr. Lisle 
Mr . Teague 
Mr . Cummings 
Dr. Sperling 
Dr. Heybey 
Mr. Mabry 
Mr. Peas ley 
Mr. 0. E. Smith 
Mr. Dickpv 


R-COMP 

Dr. Arens torf 
Mr. Davidson 
Mr. Harton 
Mr. Scollard 
Mr. Seely 
Mr. Calhoun 


22 



EXTERNAL DISTRIBUTION 


Mr. Hans K. Hinz 
Research Dept. 

Grumman Aircraft tng. Corp. 
Bethpage, Long Island, N. Y. 

Mr. J. S. Farrior 
Lockheed 
P. 0. Box 1103 
West Sta. 

Huntsville , Alabama 

Dr. George Nomicos 
Applied Mathematics Sec. 

Applied Research & Development 
Republic Aviation Corp. 
Farmingdale, Long Island, N. Y. 

Mr. Theodore E. Edelbaum 
Analytical Mechanics Assoc. Inc. 
872 Mass . Ave. 

Cambridge, Mass. 02139 


Mr. Arthur C. Gilbert, Sc. D. 
Chief, Space Systems Requirements 
Corporate Systems Center 
Div. United Aircraft Corp. 

1690 New Britian Ave. 

Farmington, Conn. 

Mr. Robert A. Lerman 
Sr. Analytical Engineer 
Tech. Planning 
Hamilton Standard Div. 

United Aircraft Corp. 

Windsor Locks, Conn. 

Dr. M. G. Boyce 
Dept, of Mathematics 
Vanderbilt University 
Nashville, Tenn. 

Dr. Daniel Dupree 
Dept, of Mathematics 
Northeast La. State College 
Monroe, La. 


Dr. Steve Hu 
Northrop Corp. 

Box 1484 
Huntsville, Ala. 

Dr. William Nesline, Jr. 

Analytical Research Dept. 

Missile and Space Div. 

Raytheon Company 
Bedford, Mass. 

Dept, of Mathematics 
Southern Illinois University 
Carbondale, 111. 

Attn: Dr. Robert W. Hunt 

Mr. Robert Silber 

Dr. H. Hermes 
Mail #A127 
Martin Company 
P. 0. Box 179 
Denver 1, Colorado 

Mr. R. W. Reck 

Martin Company 

3313 S. Memorial Parkway 

Huntsville, Alabama 

Dr. M. L. Anthony 
Space Flight Technology 
Mail No. A- 15 3 
The Martin Company 
Denver 1, Colorado 

Mr. Samuel Pines 

Analytical Mechanics Associates,* Inc. 
941 Front St. 

Uniondale, New York 

Dr. W. A. Shaw 

Mechanical Engineering Dept. 

Auburn University 
Auburn, Ala. 


23 


EXTERNAL DISTRIBUTION (Continued) 


Mr. J. W. Hanson 
Computation Center 
University of North Carolina 
Chapel Hill, North Carolina 

Mr. Oliver C. Collins 
Advanced Research 
Flight Technology Dept. 
Aero-Space Div. 

Organization 2-5762 
Mail Stop 15-12 
Boeing Co. 

P. 0. Box 3707 
Seattle, Wash. 

Mr. Carl B. Cox, Chief 
Flight Technology Dept. 
Advanced Research 
Aero-Space Div. 

Organ izat ion 2-5762 
Mail Stop 15-03 
Boeing Co. 

P. 0. Box 3707 
Seattle 24, Washington 

Mr. Richard Hardy 
Box AB-38 
Boeing Co. 

P. 0. Box 1680 
Huntsville, Ala. 

Mr. Wes Morgan 
Box AB-49 
Boeing Co. 

P. 0. Box 1680 
Huntsvil le , Alabama 

Mr. Robert Glasson 
Bendix Systems Div. 

Bendix Corp. 

3322 Memorial Parkway So. 
Huntsvil le , Alabama 


Chrysler Corp. Missile Div. 

Sixteen Mile Rd. and Van Dyke 
P. 0. Box 2628 
Detroit 31, Mich. 

Attn: Mr. T. L. Campbell 

Dept. 7162 
Applied Mathematics 

Mr. Harry Passmore 
Hayes International Corp 
P. 0. Box 2287 
Birmingham, Ala. 

As trodynamics Opera t ion 
Space Sciences Lab. 

Missile and Space Vehicle Dept. 
General Electric Company 
Valley Forge Space Technology Center 
P. 0: Box 8555 

Philadelphia 1, Pennsylvania 
Attn: Mr. J. P. deVries 

Mr. C. Cavoti 


Dr. I. E. Perlin 
Rich Computer Center 
Ga . Institute of Technology 
Atlanta, Ga . 

Dr. Charles C. Conley 
Dept, of Mathematics 
University of Wisconsin 
Mad is on, Wiscons in 

Dr. 0. R. Ainsworth 
Dept, of Mathematics 
University of Alabama 


24 



EXTERNAL DISTRIBUTION (Continued) 


Dr. N. N.Puri 

Electric Engineering Dept. 

Drexel Institute of Technology 
Philadelphia , Pennsylvania 

University of Kentucky 
College of Arts and Sciences 
Dept, of Mathematics and Astronomy 
Lexington, Kentucky 
Attn: Dr. J. C. Eaves 

Space Sciences Laboratory 
Space and Info. Systems Div. 

North American Aviation, Inc. 

12214 Lakewood Blvd. 

Downey, Calif. 

Mr. J. W. Scheuch 
P. 0. Box 557 

No. American Aviation, Inc. 
Huntsville, Ala. 

Mr. Joe Mason 

Space 6c Info. Systems Div. 

Dept. 41-595-720 

No. American Aviation, Inc. 

12214 Lakewood Blvd. 

Downey. Calif. 

Dr . D. M. Schrello, Director 

Flight Sciences 

Space 6c Info. Systems Div. 

No. American Aviation, Inc. 

12214 Lakewood Blvd. 

Downey, Calif. 

Mr. Myron Schall 

Supervisor, Powered Trajectory 

Space 6c Info. Div. 

No. American Aviation, Inc. 

12214 Lakewood Blvd. 

Downey, Calif. 


Mr. S. E. Cooper 
Space 6c Info. Div. 

Dept. 41-697-610 

No. American Aviation, Inc. 

12214 Lakewood Blvd. 

Downey, Calif. 

Space Sciences Lab. 

Space 6c Info. Systems Div. 

North American Aviation, Inc. 

12214 Lakewood Blvd. 

Downey, Cal if. 

Attn: Mr. Paul DesJardins 

Mr. Harold Bell 

Dr. E. R. Rang 
Military Products Group 
Aeronautical Div. 

Minneapolis-Honeywell Regulator Co. 
Mail Stop #340 
Minneapolis , Minnesota 

Dr . Dahlard Lukes 

M.P.G. - Aeronautical Div. 

Minneapolis-Honeywell Regulator Co. 

2600 Ridgeway Road 

Minneapolis , Minnesota 

Dr. Donald Greenspan 
Mathematics Research Center 
U.S. Army, University of Wisconsin 
Madison, Wisconsin 

Mr. Ralph Q. Haumacher 
A2-863: Space/Guidance 6c Control 

Douglas Aircraft Corporation 
3000 Ocean Park Blvd. 

Santa Monica, California 

Mr. F. A. Hewlett 
Mgr. of Documentation 
Federal Systems Div., IBM 
6702 Gulf Freeway 
Houston 17, Texas 


25 



EXTERNAL DISTRIBUTION (Continued) 


Mr. Howard S. London 
Rm. 910-E 
Bellcomm, Inc. 

1100 17th St., N. W. 

Washington, D. C. 

Mr. George Cherry 

Mass. Institute of Technology 

Cambridge, Mass. 

Mr. William C. Marshall 
Research Engineer 
M. P. G. Research Dept. 

Sta. #340 

Minneapol is-Honeywel 1 Regulator Co. 
2600 Ridgeway Roa-d 
Minneapolis 3, Minnesota 

Mr. W. G. Melbourne 
Jet Propulsion Lab. 

4800 Oak Grove Dr. 

Pasadena 3, Calif. 

Dr. Byron D. Tapley 

Dept, of Aerospace Engineering 

University of Texas 

Austin, Texas 

Dr. Hans G. Baussus 
USA, GIMRADA 
Fort Belvair, 

Virginia 


Dr. William A. Mersman 
Chief , 

Electronic Machine Computing Br. 
Ames Research Center 
Moffett Field, Calif. 

Siegfried J. Gerathewohl 
Code SB 

NASA Headquarters 
Washington, D. C. 


Dr. Herman M. Dusek 

A. C. Research & Development 

A. C. Spark Plug 

Electronics Div. of Gen. Motors 

950 N. Sepulveda Blvd . 

El Segundo, Calif. 

Mr. Howard Haglund 
Jet Propulsion Lab. 

4800 Oak Grove Dr. 

Pasadena 3, Calif. 

Mr. Hewitt Phillips 
Langley Research Center 
Hampton, Va . 

Mr. Frank J. Carroll 
Systems Requirement Dept. 
Equipment Div. 

Raytheon Co. 

40 Second Ave. 

Waltham, Mass. 

Mr. T. Perkins 
Chrysler Corp. 

HIC Bldg. 

Huntsville, Ala. 

Mr. Robert Allen 

Mgr., Huntsville Sales Office 

A. C. Spark Plug 

Electronics Div. of Gen. Motors 

Holiday Office Center 

Huntsville, Ala. 

Mr. Dale B. Ruhmel 
Staff Associate 

Aeronautical Research Assoc, of 
Princeton, Inc. 

50 Washington Rd . 

Princeton, N. J. 

Mr. Raymond Rishel 
Mathematical Analysis Staff 
Organ iza tion 2-5 330 
Mail Stop 89-75 
Boeing Co. 

P. 0. Box 3707, Seattle, Wash. 


26 



EXTERNAL DISTRIBUTION (Continued) 


Dr. Rudolf Hermann 
Director 

Univ. of Ala. Research Inst. 

4701 University Ave., N. W. 
Huntsville, Ala. 

Dr. S. H. Lehnigk 
Physical Sciences Lab. 

Army Missile Command, Bldg. 5429 
Redstone Arsenal , Ala . 

Mr. R. J. Hayes 
Code RE-TG 
NASA Hdqs . 

Washington, D. C. 

Mr. Harold Chestnut 
1 River Rd. 

Schenectady, N. Y. 

Space Sciences Lab. 

Space & Info. Systems 
No. American Aviation, Inc. 
Downey, Calif. 

Attn: Mr. Dave Engles 

Dr . Lu Ting 

Dept, of Applied Mechanics 
Polytechnic Inst, of Brooklyn 
333 Jay St. 

Brooklyn 1, N. York 

NASA, Langley Research Center 
Hampton, Va. 

Attn: Librarian 


Dr. Peter Musen 

Goddard Space Fit. Center 

NASA 

Greenbelt, Md. 

Dr. V. Szebehely 
Yale University Obs. 

Box 2023, Yale Sta. 

New Haven, Conn. 


Dr. Rudolph Kalman 

Division of Engineering Mechanics 

Stanford University 

Stanford, California 94305 

Mr. Ken Kissel 
Aeronautical Systems Div. 

Appl ied Ma thema tics Research Br . 
Wright-Pa tterson Air Force Base 
Dayton, Ohio 

Mr. Jack Funk 
Manned Spacecraft Center 
Fl ight Dynamics Branch 
NASA 

Houston, Texas 

Dr. J. B. Rosser 
Dept, of Mathematics 
Cornell Univ. 

Ithaca, N. Y. 

Douglas Aircraft Corp. 

3000 Ocean Park Blvd. 

Santa Monica, Calif. 

Attn: R. E. Holmen A2-263 

Guidance 6c Control Section 

Dr. Joseph F. Shea 

Manned Spacecraft Center (MSC) 

Houston, Texas 


Mr. Charles F. Pontious 
Guidance 6c Navigation Program 
Office of Advanced Research 6c Tech 
Code: REG 

NASA 

Washington 25, D. C. 

Mr. Jules Kanter 

Guidance 6c Navigation Program 

Office of Advanced Research 6c Tech 

Code: REG 

NASA 

Washington 25, D. C. 


27 



EXTERNAL DISTRIBUTION (Continued) 


Dr. Ray Wilson 
OART-Code RRA 
Washington 25, D. C. 

Dr. Joseph W. Siry 

Theory & Analysis Office (547) 

Da ta S ys terns D iv . 

Goddard Space Fit. Center 
Greenbelt, Md. 

Douglas Aircraft Corp. 

3000 Ocean Park Blvd. 

Santa Monica, Calif. 

Attn: Mr. Joe Santa, A2-863 

J. B. Cruz, Jr. 

Research Associate Professor 
Coordinated Science Lab. 
Urbana, 111. 

Mr. Alan L. Friedlander 
Research Engineer 
Guidance and Control Section 
IIT Research Center 
10 W. 35th St. 

Chicago 16, 111. 

Mr. Joseph V. Natrella 
Manned Space Flight - Code MCR 
NASA Hdqs . , FOB 10B 
Washington, D. C. 

Mr. Donald Jezewski 
Guidance Analysis Br. 
Spacecraft Technology Div. 
Manned Spacecraft Center 
Houston, Texas 

Mr. Robert M. Williams 
Chief, Guidance Analysis 
General Dynamics/Astronautics 
Mail Zone 513-0 
P. 0. Box 166 
San Diego 12, Calif. 


Mr. M. D. Anderson, Jr. 

Genera 1 Dynamics Corp . 

Suite 42 Holiday Office Center 
South Memorial Parkway 
Huntsville, Ala. 

Mr. Y. L. Luke 
Mathematics & Physics Div. 

Midwest Research Institute 
425 Volker Blvd. 

Kansas City 10, Missouri 

Mr. Ted Guinn 
Advanced Space Technology 
Engineering Research 
Douglas Aircraft Corp. 

Santa Monica, Calif. 

Mr. Robert Chilton 

NASA Manned Spacecraft Center - EG 

P. 0. Box 1537 

Houston, Texas 

Dr. S. E. Ross 

Office of Manned Space Fit. - NASA 
Washington, D. C. 

Lt. Stephen J. Kahne, USAF 
Applied Mathematics Branch 
Data Sciences Lab. 

Air Force Cambridge Research Labs. 
Office of Aerospace Research 
Lawrence G. Hans com Field 
Bedford, Mass. 

Mr. Daniel B. Killeen 
Computer Lab. 

Norman Mayer Bldg. 

Tulane University 
New Orleans, La. 

Dr. Bernard Friedland 
Staff Scientist - Control 
General Precision, Inc. 

Little Falls, N. J. 


28 


EXTERNAL DISTRIBUTION (Continued) 


Mr. C. H. Tross 

Mgr., Aerospace Sciences 

UNIVAC Div. of Sperry Rand Corp. 

P. 0. Box 6248, San Diego 6, Calif. 

Mr. Walter L. Portugal 
Mgr., Systems Sales 
Aerospace Group 
General Precision, Inc. 

Little Falls, New Jersey 

Mr. Ken Squires 

Goddard Space Fit. Center, 

Bldg #1, NASA 
Greenbelt, Md . 

Dr. Paul Degrarabed ian 
Astro Science Lab. 

Bldg. G 

Space Technology Lab., Inc. 

One Space Park 
Redondo Beach, Calif. 

Dr. George Leitmann 
Associate Professor, Engr. Science 
University of California 
Berkely, Calif. 

Dr. R. P. Agnew 
Dept, of Mathematics 
Cornell Univers ity 
Ithaca, N. Y. 

Dr. Jurgen Moser 

Professor of Mathematics 

Graduate School of Arts And Science 

New York University 

New York , New York 


Mr. Clint Pine 
Dept, of Mathematics 
Nor thwes tern State College 
Na tchitoches , La . 

Dr. John Gates 
Jet Propulsion Lab. 

4800 Oak Grove Drive 
Pasadena 3, Calif. 

Auburn Research Foundation 
Auburn Univers ity 
Auburn, ALa. 

Grumman Library 

Grumman Aircraft Engineering Corp. 
Bethpnge, L. I., New York 

Jet Propulsion Lab. 

Library 

4800 Oak Grove Dr. 

Pasadena 3, Calif. 

Scientific and Tech. Info. Facilit 
Attn: NASA Rep. (S-AK/RKT) 

P. 0. Box 3700 
Be thes da , Md . 

NASA Ames Research Center 
Mountain View, CaLif. 

Attn: Librarian 

NASA Flight Research Center 
Edwards Air Force' Base, Calif. 
Attn: Librarian 

NASA, Goddard Space Flight Center 
Greenbe 1 t , Md . 

Attn: Librarian 


(25) 


29 



EXTERNAL DISTRIBUTION (Continued) 


NASA Launch Operations Directorate 
Cape Kennedy, Fla. 

Attn: Librarian 

NASA Lewis Research Center 
Cleveland, Ohio 
Attn: Librarian 

NASA Manned Spacecraft Center 
Houston 1, Texas 
Attn: Librarian 

NASA Wallops Space Fit. Sta. 
Wallops Island, Va . 

Attn: Librarian 

Space Flight Library 
University of Kentucky 
Lex ing ton , Kentucky 

University of Kentucky Library 
University of Kentucky 
Lex ing ton , Ken tucky 

Office of Manned Space Flight 
NASA Hdqs . 

Federal Office Bldg. #6 
Washington 25, D. C. 

Attn: Mr. Eldon Hall 

Dr. A. J. Kelly 

Mr. Bryan F. Dorlin 
Theoretical Guidance 6c Control Br. 
NASA, Ames Research Center 
Moffett Field, Calif. 

Lt. Col. R. A. Newman 

Air Force Space Systems Div. 

M-SSVH 

Bldg. 5250 

Redstone Arsenal, Ala. 


Roger Johnson 
Honeywell Research 
Roseville, Minnesota 


30 


